Accretion column disruption in GX 1+4 
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ABSTRACT 

Daily observations of the binary X-ray pulsar GX 1+4 were made with the Rossi X- 
ray Timing Explorer (RXTE) satellite between 1997 May 16-20 as part of a four-month 
monitoring program. On May 17 the sharp dips normally observed in the lightcurve 
were all but absent, resulting in a pulse fraction / p ~ 0.5 instead of the more typical 
value of ~ 0.8 measured before and after. Also observed was a dramatic hardening of 
the 2-40 keV phase-averaged spectrum. The power-law photon index was 1.16 ± 0.02, 
whereas values of 1.6-2.0 are more typical. In terms of a Comptonization continuum 
component, the optical depth for scattering was r ~ 19, with 4-6 the usual range for 
RXTE spectra (Galloway 2000). Pulse-phase spectrosopy indicates that r is decreased 
relative to the phase-averaged value around the primary minimum, where an increase 
is normally observed. The reduced depth of the dip is interpreted as disruption of 
the accretion column, and the accompanying spectral variation suggests a substantially 
different accretion regime than is usual for this source. 

Subject headings: X-rays: stars — pulsars: individual (GX 1+4) — accretion — scat- 
tering 

1. Introduction 



The X-ray pulsar GX 1+4 is unusual in several respects. It is the only confirmed neutron 
star in a symbiotic binary (Belczynski et al. 2000) which appears to be accreting from the giant 
companions' stellar wind. The neutron star spin period evolution exhibits long periods of almost 
constant spin-up or -down, punctuated by 'torque reversals' lasting tens of days (e.g. Nelson et al. 
1997). Indirect estimates indicate a surface magnetic field strength of 2-3 xlO 13 G (Beurle et al. 
1984; Dotani et al. 1989; Mony et al. 1991; Greenhill et al. 1993; Cui 1997). These values are 
significantly greater than the ~ 10 12 G inferred for other accreting X-ray pulsars, generally from 
measurements of cyclotron scattering resonance features (CSRFs; e.g. Heindl et al. 2000) in the 
X-ray spectra. As the presence of such a strong field in GX 1+4 awaits confirmation by such 
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methods, the possibility that the field strength may in fact be one or more orders of magnitude 
lower remains. 

Historically X-ray pulsar spectra have generally been fitted with variations on a power law, 
which is well known to be a signature of unsaturated Comptonization (e.g Pozdnyakov et al. 1983). 
Model fitting to the available archival RXTE data show that the spectra of GX 1+4 are broadly 
consistent with an analytic approximation to a Comptonization continuum component combined 
with a gaussian to represent iron line emission, both attenuated by neutral absorbtion with variable 
column density (Galloway 2000). 

The archival RXTE data, which span two orders of magnitude in flux for the source, has 
also revealed that the X-ray lightcurves and mean pulse profiles of all but one of the observations 
exhibit sharp dips spanning A(f> ss 0.1-0.2 in phase. During 1996 July the countrate dropped to 
~ 10 counts -1 over ~ 6 h; sharp dips were observed in the mean pulse profiles even then (Giles 
et al. 2000). Similar dips are observed in only two other X-ray pulsars, RX J0812. 4-3114 (Reig 
& Roche 1999) and A 0535+262 (Cemeljic & Bulik 1998). A recent analytic and modelling study 
suggests that the sharp dips are signatures of eclipses of the emitting region by the accretion column 
(Galloway et al. 2000, in preparation). Pulse-phase spectroscopy using a continuum spectral model 
simulating thermal Comptonization of both GX 1+4 and RX J0812.4— 3114 shows a significant 
increase in the fitted optical depth r coincident with the dip phase. If the inclination angle i of 
these systems are roughly equal to their magnetic colatitude (3, such an increase may result when 
one of the accretion columns approaches alignment with the line of sight once each rotation period. 
At the phase of closest alignment, photons emitted from the polar cap must propagate a greater 
distance through the column to escape and ultimately reach the observer. Numerical modeling 
using a simplified column geometry confirms that the resulting increase in optical depth will have 
just the effect observed in GX 1+4 and RX J0812.4-3114 (Galloway & Wu 2000). 



2. Observations and Data Analysis 

RXTE consists of three instruments, the proportional counter array (PCA) sensitive to photons 
in the energy range 2-60 keV, the high-energy X-ray timing experiment (HEXTE) covering 16— 
250 keV, and the all-sky monitor (ASM) which spans 2-10 keV (Giles et al. 1995). Analysis of 
RXTE data presented in this letter was carried out using lheasoft 5.0, released 23 February 
2000 by the RXTE Guest Observer Facility (GOF). The data were first screened to ensure that 
the pointing offset was < 0.02° and the source was > 10° from the sun. This introduces additional 
gaps to the data. Instrumental background from cosmic ray interactions and as a result of satellite 
passages close to the SAA are estimated using the pcabackest software which is included in the 
lheasoft package. The 'bright' source background models from the 23 February 2000 release were 
used, which are appropriate when the net source countrate is > 60 counts -1 . Spectral fitting was 
undertaken using the xspec spectral fitting package version 11 (Arnaud 1996). For more details of 
the spectral analysis see Galloway (2000). 
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Between 1997 May 16-20 GX 1+4 was observed daily for ss 12 min. The phase-averaged back- 
ground subtracted countrate varied from day to day, from < 100 counts -1 (with all 5 proportional 
counter units operating) on 17 May to ~ 500 counts -1 on 18 and 20 May (Figure 1). For four 
of the five observations the sharp dips which form the primary minimum in the pulse profile were 
clearly observed in the background subtracted PCA countrate. The pulse fraction is defined as 

fp = p where 

"max 

^min < F{(t>i) < F max V i e 1 . . . n 

and where F ((/){) is the mean pulse profile calculated over n bins. On May 16 and 18-20 f p « 0.8, 
as is typical for the source. On 17 May however, the sharp dips in the pulse profile were weaker or 
absent, resulting in a pulse fraction / p « 0.5 (Figure 1, second panel from top). The minima are 
instead broad and irregular, lacking the consistency over successive cyles observed on other days. 

Spectral model fitting of the Standard-2 mode spectra averaged over each daily observation 
was undertaken with the Comptonization continuum model ('compTT' in XSPEC; Titarchuk 1994) 
following Galloway (2000). Note that it was not possible to obtain a background spectrum for 
the May 20 observation due to problems with the filter file; the background spectrum from the 
previous days' observation was used instead. Several of the spectral fit parameters show significant 
variability (Table 1). The column density «h varied between « 10 23 cm -2 and ps 3.2 x 10 23 cm -2 . 
The source spectrum temperature To was somewhat higher than is typical on 16 May but decreased 
to a more typical 1.2—1.3 keV by 18 May. For two of the spectra the scattering plasma temperature 
T e similarly could not be constrained by the fitting algorithm, and instead was fixed at 10 keV. 
For the other three days T c varied between ~ 6.5-7 and ~ 10 keV. The most dramatic variation 
was in the fitted optical depth r, which was 3.6-4 on May 16 and 18-20 but « 19 on May 17. The 
equivalent width of the Fe gaussian component was also unusually high for this spectrum, at more 
than 1.8 keV. These variations occurred as the luminosity varied between 1.5-3.6 x 10 37 ergs -1 , 
assuming a source distance of 10 kpc. The actual distance of the source is thought to be in the 
range 3-15 kpc (Chakrabarty & Roche 1997). 

The va l ues indicate good fits for each spectrum except on May 17. The spectrum obtained 
on this day does not resemble any of the other ~ 50 spectra extracted over the course of the RXTE 
observations between 1996 and 1997. The model fit was problematic in several respects. Initial 
fits with all parameters free to vary resulted in a fit value for the r max parameter (associated with 
the Xe absorption edge at 4.83 keV) an order of magnitude greater than the maximum for all the 
other spectra. Thus the value for this parameter was instead fixed ('frozen') at the fit value for the 
spectrum extracted over the whole of the 1997 May observation. The fitted To value, when free to 
vary, also became much larger than was consistent than for the other spectra and so was fixed at 
1.3 keV. This value is typical for the source at other times (Galloway 2000). While the resulting 
fit statistic xl ~ 2.60 indicates an unacceptable fit, no improvement was found using alternative 
continuum models, including power law, power law with exponential cutoff and broken powerlaw. 
A comparable fit was obtained with a blackbody component resulting in xl = 2.74 with fitted 
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T bb = 6.65|^ keV in rough agreement with T c for the compTT fit of 6.5l|;^ keV. 

The data from each of the observations on May 16-20 were divided into 10 phase bins based on 
the pulse period P = 126.4319(0) ± 0.0005(5) s estimated from regular BATSE monitoring of the 
source. The ephemeris for each day was chosen arbitrarily so that the first phase bin was centered 
on the primary minimum, defined as phase 0.0. Each of the resulting phase-selected spectra were 
fitted with the usual spectral model with at most T e , r and the compTT and gaussian component 
normalisations free to vary. The other parameters were fixed at the fit values for the phase-averaged 
spectrum on the same day. This approach is essentially identical to that of Galloway et al. (2000). 
For the May 16 and 19 observations T e was also fixed at the same value as for the phase- averaged 
spectrum (see Table 1). Variation of T c and r with phase is small except close to the primary 
minimum. The fitted r around phase 0.0 on May 16 and 18-20 was significantly greater than the 
value for the phase-averaged spectrum, by 20-60 per cent. For the May 17 observation, at the 
much shallower primary minimum r is instead marginally lower than the phase-averaged fit values, 
while T e is somewhat higher (Figure 2). 



3. Discussion 

For sources with i « /3, assuming that i, (3 and the footpoints of the accretion column on the 
neutron star do not change, dips due to column eclipses must always be observed while accretion 
continues through the columns. The apparent disappearance of the dip during the 1997 May 17 
observations clearly challenges this picture. The flux at this time is more than an order of magnitude 
larger than the lowest level measured by RXTE; complete cessation of accretion seems unlikely. 
Another possibility is that the geometry of the system is changing such that the accretion column 
drops out of alignment briefly. The inclination angle i may vary due to geodetic precession of 
the pulsar, as is suggested for some radio pulsars (e.g. Kramer 1998). Alternatively, the magnetic 
colatitude [3 may change due to evolution of the neutron star magnetic field. Variations due to the 
former effect, however, are only expected on timescales of years. And while the latter possibility 
cannot be discounted due to the lack of knowledge regarding such evolution, it seems much more 
likely that each of these mechanisms would cause a much longer-lasting, if not permanent, cessation 
of eclipses rather than what is observed in GX 1+4. 

Variations in the position of the column itself are much more plausible. As the accretion rate 
M (and hence the luminosity Lx) decreases, the balance of magnetic and gas pressures at the inner 
edge of the accretion disc means that the disk will be truncated further out from the neutron star 
(Ghosh & Lamb 1979a, b). Field lines threading the boundary region where disk material decelerates 
and becomes bound to the accretion column will thus meet the neutron star at progressively higher 
latitudes on the star. In principle at least this mechanism may cause the column to move in and 
out of alignment, in which case eclipses (and dips) would be expected only over a finite range of 
Lx- Instead what appears to be the case for GX 1+4 is that dips are normally observed at the 
full range of luminosities at which the source was observed by RXTE, and are perhaps absent at a 
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particular flux level only within that range. Observations at flux levels comparable to that of 1997 
May 17 however do generally exhibit sharp dips. If the accretion column has indeed been disrupted 
rather than moved, accretion must be occurring over a much larger region of the neutron stars' 
magnetosphere to give the observed X-ray flux. That material can enter the magnetosphere of the 
star without being channeled along the field lines suggests either that the magnetic field strength 
is in actual fact low for an X-ray pulsar, or that some global magnetohydro dynamic instability is 
operating which can completely disrupt the flow. 

The unique phase-averaged RXTE spectrum observed coincident with the unusually low pulse 
fraction underlines the significance of this event. That the optical depth has increased fivefold 
compared to values obtained before and after suggests that the source has been 'smothered' by 
hot accreting material, which presumably also contributes to the unusually strong Fe line emission. 
Such a large r suggests that the photon spectrum will approach thermalization with the scattering 
material; the close agreement of the fitted blackbody temperature and the scattering plasma tem- 
perature provides corroboration for the compTT component fits during this interval. The pulse-phase 
spectral variation was also different to that normally observed in the source. Typically r is signif- 
icantly greater than the mean value around the primary minimum; this is the case for previously 
published results (e.g. Galloway et al. 2000) as well as the observations on May 16 and 18-20. On 
May 17 r was instead somewhat less than the phase-averaged value around the primary minimum. 
These results are strongly suggestive of a substantially different distribution of matter around the 
source. The increase in the equivalent width also observed on May 17 was not accompanied by 
the usual increase in the neutral column density nu (Kotani et al. 1999; Galloway 2000), further 
indicating varying conditions close to the star rather than in the outer circumstellar matter. 

At least one similar event has been observed previously from this source. On 1993 December 
11 a pulse fraction ~ 35 per cent (20-100 keV) was measured using balloon-borne instruments, 
along with an unusually hard spectrum with power law energy index 0.54 ± 0.18 (Rao et al. 1994). 
For comparison, a power law fit was made to RXTE spectra from 17 May in the energy band 
20-40 keV only. The resulting photon index 1.15} o| is consistent with the fit value for the full 
PC A energy range, and represents a significantly harder spectrum than on 1993 December 11. 
Flux measurements by the Burst and Transient Source Experiment (BATSE) aboard the Compton 
Gamma-ray Observatory satellite indicate that following that observation the source entered a hard 
X-ray low state with Lx < 0.1 keVcm~ 2 s (20-60 keV) which persisted (with brief interruptions) 
for the next s» 250 d. In contrast, the May 17 observation was made during a flare lasting ~ 20 d 
where the flux peaked at 0.2 keVcm~ 2 s . The other spectra from 1997 May give fitted r values 
at intermediate levels between the two distinct spectral states suggested by analysis of the archival 
RXTE observations (Galloway 2000). It is possible that the disruption of the accretion column is 
also related to transitions between distinct spectral and luminosity states inferred from analysis of 
the full set of RXTE archival data (Galloway 2000). 

This research has made use of data obtained through the High Energy Astrophysics Science 
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Archive Research Center Online Service, provided by the NASA/Goddard Space Flight Center, 
and also the BATSE pulsar group home page at http://www.batse.msfc.masa.gov. 
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Fig. 1. — RXTE observations of GX 1+4, 1997 May 16-20. Each panel shows the countrate on 1 s 
bins for all 5 PCUs for the short observations on May 16, 17, 18, 19 and 20 from top to bottom 
respectively. Shown at the right are the corresponding mean pulse profiles (over two full cycles) 
with a representative error bar. The profiles are calculated by folding the 1 s lightcurves on the 
pulse period 126.4319(0) ± 0.0005(5) s estimated from regular BATSE monitoring of the source. 
Note the change in scale on the y-axis between the panels. 
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Fig. 2. — Pulse-phase spectroscopy of data from the RXTE observation of GX 1+4 on 1997 May 17. 
The top and middle panels show the variation in fitted T e and r respectively with phase. The solid 
and dotted lines show the fitted values and 90 per cent uncertainty limits for the phase-averaged 
spectrum. Errorbars represent 90 per cent uncertainty limits. The bottom panel shows the pulse 
profile over the full PCA energy range (as shown in Figure 1). 
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Table 1. Spectral fit parameters using the Comptonization continuum component for the May 
1997 RXTE observations of GX 1+4. Fit parameters are the absorption column density nn (in 
units of 10 22 cm -2 ), source spectrum temperature To, plasma temperature T e and optical depth 
r, Comptonization component normalisation Aq (in units of 10~ 2 photons cm~ 2 s -1 keV -1 ), 
gaussian component center energy E Fe , width a, normalisation ^4f c (in units of 
10~ 3 photons cm" 2 s -1 keV" 1 ) and equivalent width EW, and finally the 2-60 keV flux Lx in units 
of 10 37 ergs -1 (assuming a distance of 10 kpc). Energy range for fitting is typically 2.8-35 keV. 



Parameter 
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